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ABSTRACT 



Context. AADor is an eclipsing, post common-envelope binary with an sdOB-type primary and a low-mass secondary. Eleven years 
ago, an NLTE spectral analysis showed a discrepancy in the surface gravity that was derived by radial-velocity and light-curve 
analysis, log^ = 5.21 ± 0.1 (cm/sec-) and logg = 5.53 ± 0.03, respectively. 

Aims. We aim to determine both the effective temperature and surface gravity of AADor precisely from high-resolution, high-S/N 
observations taken during the occultation of the secondary. 

Methods. We calculated an extended grid of metal-line blanketed, state-of-the-art, non-LTE model atmospheres in the parameter 
range of the primary of AADor. Synthetic spectra calculated from this grid were compared to optical observations. 
Results. We verify Teff = 42000 ± lOOOK from our former analyses and determine a higher log^ = 5.46 ± 0.05. The main reason are 
new Stark-broadening tables that were used for calculating of the theoretical Balmer-line profiles. 

Conclusions. Our result for the surface gravity agrees with the value from light-curve analysis within the error limits, thereby solving 
the so-called gravity problem in AA Dor. 
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1. Introduction 

AA Dor is a close, eclipsing, post common-envelope binary sys- 
tem with an sdOB-type primary star and an unseen low-mass 
companion. The orbital period is 0.261 539 7363 (4) d (Kilkenny 
201 1[) and the inclination is / = 89°21 + 0°30 (Hilditch et al. 
2003h . A detail ed introductio n to the system and previous anal- 
yses is given in lRauchI (l2004t) . and we summarize results of pre- 
vious spectral analyses of the primary in Table[T] 

[Rauch (2000) encountered the problem of his spectroscop- 
ically determined surface gravity logg = 5.21 +0.1 not match- 
ing log g = 5.53 ± . 03 determined from light-curve analysis 
dHilditch et al.lll996h . lHilditch et all (12003) present an improved 
photometric model and derive log;? = 5.45 - 5.51 . The reason 
for the log g discrepancy is unknown. iFleig et aH (j2008) find a 
slightly higher lo g ^ = 5.3 + 0.1 but th e discrepancy remains. A 
recent analysis by Mii ller et al.l (1201 Oh has apparently solved the 
log ;? problem by finding log ^ = 5.51 +0.05. 

iMiiller etal.1 (2010) do not consider that the He i lines 
(as well as other lines of low-ionized species, e.g. of Mg ii, 
Fig.[T]i are too strong in the models at their favored parameters 
Teff = 37800 + 500 K and log g = 5.5 1 + 0.05, as demonstrated in 
FigHI Increasing the He abundance to better fit He ii /1 4686 A 
also results in a much stronger He i A 447 1 A, which then dis- 
agrees with the observation. We mention that Fleig et all (l2008l) 
evaluate the ionization equilibria of C m / C iv, N iii / N iv, O iii 
/ O IV, P IV / P V, and S iv / S v in the FUV wavelength range 
and find Tetj = 42000 ± 1000 K to agree with the higher Teff con- 
cluded from the He i lines. 
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Fig. 1. Comparison of our synthetic spectra (full, blue 
line: reff = 37800K, logg = 5.51; dashed, red: refF = 42000K, 
log § = 5.46; He = 0.0027 by mass) around He 1/14471 A and 
Mg II /1 4481 A (left) and He ii /1 4686 A (right) with the obser- 
vation. The models are convolved with a rotational profile cor- 
responding to v,ot = 30 km/sec. Models and observation are 
smoothed with a Gaussian (0. 1 A FWHM) for clarity. 



* Based on observations made with ESO Telescopes at the Paranal 
Observatories under programme ID 66.D-1800. 



Since iKuruczl (^2009', http://kurucz.harvard.edu/atoms.htm^ 
has substantially extended his database, and the model atoms in 
our Tubingen Model-Atom Database (TMAZfl) have been up- 
dated as well, we decided to calculate an improved, extended, 
state-of-the-art NLTE model-atmosphere grid. This grid is de- 
scribed in Sect.|2] The re-analysis of our UVES spectra (105 a 
180 sec, which in total cover one orbital period and which were 



http://astro.uni-tuebingen.de/~TMAD/TMAD.html 
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Table 1. Effective temperature and surface gravity of the primary of AADor, determined in previous and the present spectral 
analyses. 



^eff 




He 


I'm! 






method 


reference 


(K) 


(cm/sec^) 


(mass fraction) 


(km/sec) 


{Me) 


(Me) 






41000 


5.4 


0.28 








LTE' 


Kudritzki (1976) 


44200 


5.2 


0.28 








NLTE' 


Kudritzki (1976) 


41700 


5.9 


0.8 








LTE' 


Kudritzki (1976) 


42000 


5.7 


0.8 








NLTE' 


Kudritzki (1976) 


40000+'"'™ 

tuuuu 2000 


5.3 ±0.2 


0.012 




0.3 


0.04 


NLTE' 


Kudritzki et al (1982) 


42000±1000 


5.21±0.1 


0.0032 


34.0 


0.33 


0.066 


NLTE^ 


R.auch (2000)" 


42000±1000 


5.30±0.1 


0.0032 


35.0 






NLTE^ 


FleiH et al. (2008) 


37800± 500 


5.51±0.05 


0.005 


30.0 


0.51 


0.085 


LTE'' 


Muller etal. (2010) 


42000±1000 


5.46±0.05 


0.0027 


30.0 


0.47 


0.079 


NLTE' 


this work 



' H+He models, two grids with fixed NaJNii ratios only, 1"' investigation of NLTE effects, no errors given 

^ H+He+C+N+O+Mg+Si+Fe+Ni, Fe+Ni data from lKuruc3 <1991l) . optical spectra, assumed bound rotation (vf^^ = 45.7 km/sec) 
3 H+He+C+N+O+Mg+Si+P+S+Ca+Sc+Ti+V+Cr+Mn+Fe+Co+Ni, Ca-Ni data fromJCurucz ( 1991), optical and FUV spectra 

H+He metal enhanced (z = ten times solar), optical spectra 
' H+He+C+N+O+Mg+Si+P+S+Ca+Sc+Ti+V+Cr+Mn+Fe+Co+Ni, Ca-Ni data from lKurucd <2009h . optical spectra, vf^| from lMuUer et al.l ilQld) 



also used by iMiilleret all 1201 Oh that were obtained in 2001 at 
the VLT is described in Sect.|3] We conclude in Sect.|4] 



2. Atomic data and model-atmosphere grid 

The model atmospheres used here were calculated w i th the 
Tubingen Model-Atmosphere Package dWerner et al.l l2003l 
TMAP). The models are plane-parallel, in hydrostatic and radia- 
tive e quilibrium. TMAP uses the occupation-probability formal- 
ism oflHummer & Mihalas (1988) that was generalized to NLTE 
conditions by Huben y et al. (1994^). TMAP considers opaci- 
ties of HH-HeH-C-H N-nO+Mg+SiH-P-i-S using classical model 
atoms, and Ca-nSc+li+V+Cr+Mn-i-Fe-i-Co-i-Ni uses a statisti- 
cal approach (Rau ch & Dee tien 2003). All model atoms used in 
our calculations were updated to the most recent atomic data 
(Sect.[T]i, and 530 levels are treated in NLTE with 771 individual 
lines (fror n H - S) and 1 9 957 605 lines of Ca-Ni from Kurucz' 
line lists (iKurucj 12009 ) combined to 636 superlines. The ele- 
ment abundances are summarized in Table|2] 

The model-atmosphere grid spans Tefj - 35000 - 49000 K 
(ArefF = 500K) and logg = 5.15 - 6.20 (Alog§ = 0.05). In total 
this makes 638 models. Spectral energy distributions (SEDs) 
were calculated using the most recent line broad ening data, 
e.g. H I line-broadening has changed in TMAP since iFleig et all 
(l2008h presented the ir analysis of AADor. The reason is that 
iRepolust et al.l (12005*) found an eiTor in the H i line-broadening 
tables (for high members of the spectral series only) by Lemkg 
(Il997h that were used before. These were substituted by a 
Holtsm ark approximation. In addition, Tremblay & Berg erog 
( |2009|) provide new, parameter-free Stark line-broadening tables 
for Hi considering non-ideal effects. These replaced Lemke's 
data for the lowest ten members of the Hi Lyman and Balmer 
series. In the parameter range of AADor, the new broadening 
tables have a significant impact on the line wings of higher 
Balmer-series members (narrower for H e and higher, Fig.|2]i. As 
a consequence, our analysis results in a higher log g (Sect.[3]l. 

In the framework of the Virtual ObservatorjQ {VO), all these 
SEDs {A - F^) are available in VO compliant form via the VO 




4000 4500 

1/A 

Fig. 2. Synthetic spectrum calculated from a T^ff - 42000 K and 
log g = 5. 45 mode l with d ifferent Stark line-broadening tables (L , 
blue line: iLemkd (|1997|) . T, red: iTremblav & BergeronI (|2009|) . 

see text). 



service 77!eo5&4Q provided by the German Astrophysical Virtual 
Observatory (GAVCQ). 



3. Analysis and results 

The light curve of AADor exhibits a reflection effect (e.g. 
Hilditch et al. 1996) that amounts to about 0.06 mag in the op- 
tical. To analyze the pure primary spectrum, we selected only 
those four observations that were taken closest to the occulta- 
tion of the secondary. These were co-added in order to improve 
the S/N. In Fig.|3] we show a fit to all single UVES spectra. 
Our fit excludes the inner line cores of H/? and Hy, as well 
as obviously bad data points (quality flags given bv lMiiller et al.l 
2010). Both the occultation (at (j) - 0.5) and the transit (0 = 0.0) 
of the secondary are clearly visible in the determination of T^g 
and logg. Compared to a similar x^ fit of iMiiller et alJ (1201 Ol 



^ http://www.ivoa.net 



^ http://vo.ari.uni-heidelberg.de/ssatr-0.01/TrSpectra.jsp? 
* http://www.g-vo.org 
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Fig. 4. Left: Formal 1 cr, 2 cr, and 3 cr contour lines of our fits 
in the T^s - log g plane. Right: Reduced of our T^s - 42000 K 
models depending on log g. 



their Fig. 3), we find the same log§ = 5.45 but a significantly 
higher Teft =40500 K than for T^s = 37800 K. 

For the analysis, we perform a detailed comparison in the 
classical way (;t^-b y-eye) and, for comparison in analogy to 
iMtiller et alJ(l2010l) " with ?lx fit, used the same wavelength lim- 
its (Tables and lines, H/3 - Hll and Hen ^4686 A. Om 
fit yields Teff = 40600 + 100 K and log g = 5.46+J}j;^ (T in Fig.©. 
These errors are formal 1 cr errors, and cr was calculated from 
the deviation of thej^f^^j^^ model from the observed spectrum used 
in the x^ fit- Compared to a similar x^ fit with SEDs that were 
calculated with the previously used Stark broadening tables of 
iLemke, (.1997. L in Fig.|4]i, there is a significant deviation of 
ATeff = 600 K and A log g = 0.06. 

A comparison of the best-fitting model from our x^ fit and 
the best-fitting ;(f-by-eye with the observations is shown in 
Fig. El It is obvious that the ionization equilibrium of He i/He ii 
is reproduced not at reff = 40600 K but at reff = 42000 K. The 
theoretical line profiles of lower members of the Balmer se- 
ries (Hy6 - H(J) do not reproduce the observation perfectly. 
They fit sHght ly better at r^ff^ 40600 K. T hus, a small Balmer- 
line problem ( iNapiwotzki & Rauc"hlll994l: IWern er 1996) due to 
additional metal opacities that are still not considered is ap- 



0.5 



0) 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

Fig. 3. Phase-dependent, best-fitting grid model determined by a 0.0 
X^ fit. <p - 0.0 is the transit, (p - 0.5 is the occultation of the sec- 
ondary. (The steps in Tetf and \ogg represent the grid spacing.) 
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Fig. 5. Comparison of synthetic line profiles of H and He lines 
calculated from a Teff^ 40600 K and log g = 5.46 model (left) 
and a Teff - 42000 K and log § = 5.46 model (right) with the ob- 
servation. 



parently present. The inclusion of He 1/14471 A (TableO in 
the x^-^^ procedure results in higher re ff= 40700 + 300 K . We 
finally adopt Teff =42000 + 1000 K (cf. Fleig et al. 2008) and 
log ^ = 5.4 6 ± 0.05 because the previous ly evaluated ionization 
equilibria (lRauchll2000HFleig et al.ll2008h are an additional, cru- 
cial constraint. Kx^ fit at fixed T^ti - 42000 K (additional models 
were calculated with logg = 5.30 - 5.60 and Alog§ = 0.01) also 
has its minimum at log g = 5.46 (Fig.|4li. 



A mass of M, 



pri 



0.4714 + 0.0050 Mo is determined by 
comparing of Tetf and log g with the evolutionary tracks of post- 
EHB stars (Fig.|6]l. From the same evolutionary calculations, we 
interpolate the primary's luminosity. From our final model, we 
can determine the spectroscopic distance of AA Dor following 
iHeber at al.1 (11984 . We derive a distance of c/ = 352+?? pc. The 
parameters of AADor are summarized in Tables|2]and]4] 



Tab le 4. Parameter s of AADor compared with values of 
.Hilditch et al.l(l200l . 





this work 




Hilditch et al. f2003) 


l0g(gp,i/f ) 


42.0 
5.46 


±1 

±0.05 


5.45 -5.51 


Mp,,/Mo 


0.4714 


±0.0050 


0.33 -0.47 


^pri/ L'Q 

rf/pc 


120 
0.0788 
352 


+ 15 
-20 

+0.0075 
-0.0063 
+20 
-23 


0.064 - 0.082 



4 
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Fig. 6. Location of AADor in the Tgff-logg plane compared to 
sdBs and sdOBs from Edelmann (2003J. Post-EHB tracks from 
iDorman et alj (Il998l labeled with the respective stellar masses 
in Mq) are also shown. 



4. Conclusions 

The so-called \ogg problem in AADor is solved (Fig.|2l) and 
our resultsJTablelDi are in good agreement with the photometric 
model of imiditch et all (120031) . 

Four influences were identified on the logg determination. 
1) The major impact is the improvement in the Stark broa dening 
tables, i.e. the difference between those of Lemke ('1997) and of 
[Tremblav & Bergeron (2009). This results in a systematic devia- 
tion of AFeff = 600 K and A log ^ = 0.06. 2) The reflection effect 
is now eliminated by using only observed spectra that were ob- 
tained during the occultation of the secondary (Sect. [3]). 3) The 
improved atomic data makes the model -atmosphere more reli- 
able thanks to a fuller consideration of the metal-line blanketing. 
The temperature stratification of the stellar models, however, is 
only marginally affe cted. 4) The rotation al velocity is lower than 
previously assumed dMiiller et al.ll20ldl) . This only has little in- 
fluence on the inner line core and is thus important for weak and 
narrow li nes like He ii /1 4686 A (Fig.[I]l. 

Since IVuckovic et alj (l2008l) identified spectral lines of the 
secondary in the UVES spectra and determined a lower limit 
(/Tsec > 230 km/sec) of its orbital velocity amplitude, both com- 
ponents' masses are k nown (MpH - 0.45 Mq, Msec = 0.076 Mp, 
[Vuckovic et al."2008'), albeit with large error bars. i Miiller et alJ 
(2010) used the velocity amplitudes of both components (K^n - 
40.15 ± 0.11 km/sec, K^ec = 240 + 20 km/sec) to derive the 



masses M, 



pn 



0.51+;;;;,-;= Mo and M,^ 



0.085t°°llMQ. This 



rules out a post-RGB scenario because post-RGB masses are sig- 
nificantly lower. The solution from mass function f(m) and light 
curve analysis, however, intersects with our result of logg - 
5.46 + 0.05 (Fig.|7]) even for the higher post-EHB mass (Fig.|6l). 



From our mass determination of M, 



pn 



0.4714 



0.0050 Mq, we calculated (MpnKpn = Msec^sec) the secondary's 
mass of Msec = 0.0725 - 0.0863 Mq. Since the hydrogen- 
burning mass limit i s about 0.075 Mq (Chabrier & Baraffe 1997; 
IChabrier et al.l2000l) . the secondary may either be a brown dwarf 
or a late M dwarf. 



Fig. 7. Mass-radius relation for the primary of AADor. The 
dashed lines show the error ranges. The vertical lines sh ow the 
primary mass, derived from comparison with post-RGB (iRauchl 
i2000l) and post-EHB (Fig.|6]l evolutionary models. 
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Table 2. Element abundances in our model-atmosphere grid. 



element 


mass 


number 






X 


fraction 


fraction 




[X]" 



H 


9.939E- 


-1 


9.992E- 


-1 


12.144 


+0.130 


He 


2.686E- 


-3 


6.801E- 


-4 


8.977 


-1.967 


C 


1.777E- 


-5 


1.499E- 


-6 


6.320 


-2.124 


N 


4.145E- 


-5 


2.998E- 


-6 


6.621 


-1.223 





1.009E- 


-3 


6.393E- 


-5 


7.950 


-0.754 


Mg 


4.078E- 


-4 


1.700E- 


-5 


7.375 


-0.240 


Si 


3.049E- 


-4 


I.IOOE- 


-5 


7.186 


-0.339 


P 


5.197E- 


-6 


1.700E- 


-7 


5.375 


-0.050 


S 


3.241E- 


-6 


1.024E- 


-7 


5.155 


-1.980 


Ca 


5.994E- 


-5 


1.515E- 


-6 


6.325 


-0.030 


Sc 


3.694E- 


-8 


8.327E- 


-0 


3.065 


-0.100 


Ti 


2.784E- 


-6 


5.894E- 


-8 


4.915 


-0.050 


V 


3.731E- 


-7 


7.421E- 


-9 


4.015 


+0.070 


Cr 


1.663E- 


-5 


3.240E- 


-7 


5.655 


+0.001 


Mn 


9.877E- 


-6 


1.828E- 


-7 


5.405 


-0.040 


Fe 


1.153E- 


-3 


2.091E- 


-5 


7.465 


-0.050 


Co 


3.591E- 


-6 


6.174E- 


-8 


4.935 


-0.069 


Ni 


3.482E- 


-4 


6.009E- 


-6 


6.923 


+0.689 



log(e/6^). XoaTr.u.e, = 12.15. fcf. lHoIwegeilll979t) 
**: log[abundance/solar abundance] (solar values from lAspIund et aDl2009l) 

Table 3. Lines and wavelength intervals used for our;^^^ fits. 



line 




line 






(A) 




(A) 


H/J 


[-50,+50] 


He I ^4471 A 


[ -5, +5] 


Hy 


[-40, +40] 


He Hi 4686 A 


[ -5, +5] 


H(5 


[-30,+30] 






He 


[-20,+20] 






H8 


[-20,+20] 






H9 


[-10,+10] 






HIO 


[-10, +7] 






Hll 


[-10, +8] 







